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Abstract 

A detailed study of inclusive deep inelastic scattering {DIS) from mirror ^ = 3 nuclei 
at large values of the Bjorken variable x is presented. The main purpose is to estimate 
the theoretical uncertainties on the extraction of the neutron DIS structure function 
from such nuclear measurements. On one hand, within models in which no modification 
of the bound nucleon structure functions is taken into account, we have investigated 
the possible uncertainties arising from: i) charge symmetry breaking terms in the 
nucleon-nucleon interaction, ii) finite eff^ects neglected in the Bjorken limit, iii) the 
role of different prescriptions for the nucleon Spectral Function normalization providing 
baryon number conservation, and iv) the differences between the virtual nucleon and 
light cone formalisms. Although these effects have been not yet considered in existing 
analyses, our conclusion is that all these effects cancel at the level of ~ 1% for x < 0.75 
in overall agreement with previous findings. On the other hand we have considered 
several models in which the modification of the bound nucleon structure functions is 
accounted for to describe the EMC effect in DIS scattering from nuclei. It turns out 
that within these models the cancellation of nuclear effects is expected to occur only 
at a level of ~ 3%, leading to an accuracy of ~ 12% in the extraction of the neutron to 
proton structure function ratio at x ~ 0.7 0.8. Another consequence of considering 
a broad range of models of the EMC effect is that the previously suggested iteration 
procedure does not improve the accuracy of the extraction of the neutron to proton 
structure function ratio. 
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1 Introduction 



The investigation of deep inelastic scattering (DIS) of leptons off the nucleon is an important 
tool to get fundamental information on the structure of quark distributions in the nucleon. 
In the past years several experiments have been performed in order to study the region of 
small values of the Bjorken variable x = Q^/2Mz/, which is dominated by sea quarks and 
gluons. Recently, experiments at HERA have pushed the measurements at large x to a new 
high-Q^ domain, while dijet measurements at Tevatron have reached the kinematics where 
the knowledge of the quark distributions in the nucleon at x > 0.6 becomes important (for 
a recent discussion and references see 

One of the major uncertainty in large-x studies comes from a poor knowledge of the 
(i-quark distribution in the nucleon. The reason is that it is very difficult to extract d{x, Q"^) 
from measurements off hydrogen targets: it enters as a correction in case of inclusive electron 
scattering off the proton at <^ (where Mw is the mass of the VT-boson), while the 
measurements using semi-inclusive e + p^e + ir^ + X and large-Q^ + p —* z/(z/) + X 
scatterings do not have reached a sufficient degree of statistical accuracy yet. As a result, 
one has to rely on the extraction of d{x, Q"^) from the data involving DIS off the deuteron. 

On the theoretical side the predictions for the behavior of the ratio d{x, Q'^)/u{x, Q"^) at 
a; — > 1 vary very significantly. The deviations from the SU{6) value d{x, Q'^)/u{x, Q"^) = 0.5 
could come either from non-perturbative effects, which could lead, as suggested firstly by 
Feynman g, to d{x,Q^)/u{x,Q^) (corresponding to F:^{x,Q'^)/F^{x,Q^) 1/4) at 
s ^ 1 or from the hard scattering (perturbative QCD) mechanism, yielding d{x, Q'^)/u{x, 
1/5 (corresponding to {x , Q'^) / F2 {x , Q"^) — > 3/7 ~ 0.43) at a; — *• 1 0. Thus experimental 
data on F2{x, Q'^)/F2{x, Q^) at large x have a high degree of theoretical significance. 

The first extraction of the ratio F^{x,Q'^)/Fi{x,Q^) from the SLAC DIS p{e,e')X 
and D(e,e')X data (see [|](a) and references therein) used the West procedure [0], which is 
based on a covariant electron-deuteron scattering formalism with the interacting nucleon off- 
mass-shell and the spectator nucleon on- mass-shell. Such a procedure leads to the so-called 
West correction atotin* + d) / {atotin* +P) + ^toti.1* + t^)) < 1 in the impulse approximation 
(the numerical value of this ratio is around 0.980 -r- 0.985). The application of the West 
procedure has provided a limiting value {F2 / F2)\x-*i ~ 1/4, which has been adopted in 
most of the global fits of parton distribution functions (PDF's) (see, e.g., ||, |^, |1^, |ll[). 
Later [0, Q it was pointed out that the West correction leads to a violation of the Gross- 
Llewellyn Smith (baryon charge) sum rule, because of the neglect of relativistic corrections 
in the normalization of the deuteron wave function in Ref. 0. Furthermore it was pointed 
out in |r2| that modeling the deuteron wave function with one on-mass-shell- and one off- 
mass-shell nucleon without taking into account other degrees of freedom unavoidably leads 
to the violation of the energy-momentum sum rule, which expresses the requirement that 
the sum of the light cone fractions carried by all partons adds up to 1. Then an alternative 
light-cone (LC) formalism was suggested which satisfies both the baryon charge and 
the momentum sum rules. Both the Virtual Nucleon Convolution (VNC) model [1^ and 
the LC formalism lead to an enhancement of F^lxjQ"^) at large x as compared to the 
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SLAC procedure and hence to a further decrease of the extracted value of the n/p ratio, 
Fl'/Ff, at large x. 

The discovery of the EMC effect [15| has clearly indicated gross deviations of the jF^ 
ratio from the predictions based on the Fermi motion approximation in the kinematical 
regions at large x relevant for the extraction of F^jFl^ from the deuteron inclusive data. 
This led immediately to the conclusion |[16| that the value of F^/Ff extracted by SLAC 



was underestimated. Therefore an approximate procedure, which was argued to depend only 
marginally on the details of the EMC effect |]16|, O, was suggested for x < 0.7. The 



use of this procedure by the SLAC experimental group ||I9[ has confirmed the conclusion 
of Ref. [|l^], finding that the value of F2 /F2 at a; ~ 0.7 may be much closer to the pQCD 
asymptotic value of 3/7 ~ 0.43. 

Over the years a number of further studies of deuteron structure functions have been 
performed using the VNC model and most of them have adopted the normalization 
of the deuteron wave function of Ref. |]T2[. Some of these studies have included also pion 
degrees of freedom to fix the momentum sum rule problem. Other studies have considered 
off-mass-shell effects in the structure function of the interacting virtual nucleon. Within the 

where the 
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latter the issue of the extraction of the ratio F^/Ff' was analyzed in Ref. 
extracted values for the n/p ratio at x < 0.7 turned out to be very similar to the findings of 
Refs. P, g. 

It should be also mentioned that quite recent analyses ^ of leading and higher 
twists in proton and deuteron DIS data have found that the latter are consistent with a 
significant enhancement of the d-quark distribution at large x with respect to the standard 
PDF behavior oi d/u ^ 0. 

The realization that the extraction of the large-x n/p ratio from deuteron DIS data is 
inherently model dependent has led to the suggestion of two new strategies. One is the use 
of the tagged semi- inclusive processes off the deuteron |]23|, |2^, which require the detection 
of a low momentum proton {p < 150MeV/c). The other one is the determination of the 
DIS structure functions of mirror A = 3 nuclei In the former one can tag the 

momentum of the struck neutron by detecting the slow recoiling proton; in this way it is 
possible to select initial deuteron configurations in which the two nucleons are far apart, so 
that the struck nucleon can be considered as free. In principle, one can use here an analog of 
the Chew-Low procedure for the study of scattering off a pion and extrapolate the cross 
section to the neutron pole. The neutron structure function can then be extracted directly 
from the semi-inclusive deuteron cross section without significant nuclear model dependence 
In the latter one tries to exploit the mirror symmetry of A 



3 nuclei; in other 

words, thanks to charge symmetry, one expects that the magnitude of the EMC effect in 
'^He and 



Fi{x,Q') + F-{x,Q') 
ZFi{x,Q^) + NF^{x,Q^' 
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is very similar and hence the so called super-ratio 



should be very close to unity regardless of the size of the EMC ratios itself 26|. If 
this is true, the n/p ratio could be extracted directly from the ratio of the measurements of 
the ^He to DIS structure functions without significant nuclear modifications. However, 
it should be pointed out that, even if charge symmetry were exact, the motion of protons 
and neutrons in a non isosinglet nucleus (let's say ^He) is somewhat different due to the 
spin-flavor dependence of the nuclear force. 

The aim of this paper is to perform explicit calculations of the EMC effect for both 
^He and targets, taking properly into account the motion of protons and neutrons in 
mirror A = 3 nuclei. We explore in greater details the VNC model used in [^, ^ in 
order to analyze the effects of: i) charge symmetry breaking terms in the nucleon-nucleon 
(NN) interaction; ii) finite effects in the impulse approximation; iii) the role of 
different prescriptions for the nucleon Spectral Function normalization providing baryon 
number conservation; and iv) the role of different PDF sets. Additionally we compare 
the predictions of the VNC model and the LC formalism in the approximation where no 
bound nucleon modification is taken into account. It will be shown that the inclusion of 
these additional effects leaves the super-ratio (H) close to unity within 1% only for x < 0.75, 
confirming therefore the findings of Refs. ||25| and where deviations of the order of 2% 
and 1% were found, respectively. 

However, it is well known that the VNC model underestimates significantly the EMC 
effect at large x. Also, if the VNC model is adjusted to satisfy the momentum sum rule by 
adding pionic degrees of freedom, it leads to a significant enhancement of the qa/qn ratio at 
X > 0.1, where a suppression is observed experimentally [^. Moreover, the VNC model is 
just one of the many models of the EMC effect. Similarly the LC formalism without bound 
nucleon modifications strongly disagrees with data at large x. Hence, to provide a more 
conservative estimate of the possible range of deviations of the super-ratio from unity we 
will also investigate carefully various models of the EMC effect which interpret this effect 
as due to modification of the wave function of either individual nucleons or two nucleon 
correlations. We will show that the cancellation of the nuclear effects in the super-ratio (§) 
within the broad range of the models considered occurs only at the level of ~ 3%, restricting 
significantly (up to ~ 12%) the accuracy of the extraction of the free n/p ratio from the 
ratio of the measurements of the ^He to "^H DIS structure functions. 

In this work we will not address all the EMC models predicting possible deviations from 
the convolution formula at large x. We feel however that it's worthwhile to mention at least 
few of them. An important issue in modeling the EMC effect is the possible role of final 
state interaction effects even in the Bjorken limit. Though these effects are absent if the 
scattering process is formulated directly in terms of parton degrees of freedom, the final 
state interactions may be present in the case of a two stage descriptions, where the nucleus 
is described as a system of hadrons and next the scattering off the parton constituents of the 
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hadrons is considered. However, it is very difficult to obtain safe estimates of such an effect, 
and therefore we have not included it in the present work. Another issue is the polarizability 
of the nucleon into A(1232) components, which contributes to three-nucleon forces adding 
a ~ 10% correction to the binding energy of the three-nucleon system. Consequently, it 
may be possible that at large x one is not measuring only the nucleon structure function. 
Note here that the interference among the scattering off A-isobar and nucleon is known to 
be relevant for the description of the polarized A = 3 structure functions ||29|. The role of 
the A(1232) component effects in the problem of extraction of the F2/F2 ratio deserves a 
special study, which is beyond the scope of this paper. 

The plan of this paper is as follows. In Section 2 the formalism needed to evaluate the 
nuclear structure functions F<f^{x, Q"^) within the VNC model and the LC approximation is 
presented. All the necessary inputs for a realistic estimate of the basic nuclear ingredient, 
namely the invariant nucleon Spectral Function and its proper normalization, are discussed. 
In Section 3 both the nuclear EMC effect and the super-ratio in mirror A = ?> nuclei are 
evaluated adopting the VNC model and the LC approximation, assuming also no modifica- 
tion of the bound nucleon structure functions. Section 4 is devoted to estimate the deviations 
of the super-ratio STZemc{Xi Q"^) from unity in several models of the EMC effect, in which 
modifications of the bound nucleon structure functions are considered. The issue of the 
extraction of the n/p ratio from the measurement of the ratio of the mirror A = 3 struc- 
ture functions is fully analyzed in Section 5. Our main conclusions are then summarized in 
Section 6. 



2 Basic Formalism for Inclusive DIS from Nuclei 

There exist a number of treatments in the literature. However some of them do large 
approximations right away, do not specify completely a prescription for treating off-mass- 
shell effects in the amplitude of virtual photon-nucleon interaction, etc. Hence we find it 
necessary in this Section to rederive the basic formulae needed for the evaluation of the 
nuclear structure function F,f{x, Q"^) within the VNC model and the LC approach at finite 

In both cases no modification of the bound nucleon structure functions will be considered. 
We will refer to these approximations as convolution approximations. 

2.1 Virtual Nucleon Convolution Model 

The cross section for the inclusive A{e, e')X reaction can be written in the following general 
form 

d(J EL 

Vf^u , (3) 



dE'jn'^ Ee g4 

where rj^^ = |Tr(^27/i^i7!/) is the leptonic tensor, ki = {Ee, ki) and k2 = {E'^, k2) are the four- 
momenta of the incident and scattered electrons, respectively, and is the electromagnetic 
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tensor of the target, viz. 



wT= E {A\rA<i)\x){x\j^^{^)\A). 

spin,X 



(4) 



Within the covariant impulse approximation one assumes that the virtual photon inter- 
acts with a virtual nucleon and the final hadronic state X consists of the product of inelastic 
7*A^ interaction and the recoil [A — l)-nucleon system. Based on the Feynman diagram 



analysis of this scattering for the electromagnetic nuclear tensor ||30[ one obtains: 

Here the invariant nucleon Spectral Function in the nucleus is defined as: 

Sn{p) = / d[pA-i] T^ip, [pA~i]), 



(5) 



(6) 

where p is the momentum of the virtual nucleon, denotes internal variables of the resid- 

ual on-mass-shell {A — l)-nucleon system and r{p,pji_i) is the covariant A ^ N,{A — 1) 
vertex function combined with the propagator of the virtual nucleon. Based on the require- 
ment of baryonic number conservation the nucleon Spectral Function is normalized as 
follows: 



>4 APO 

d P -Ti— 

Ma 



Sn{p) = 1. 



(7) 



To proceed further, we express the electromagnetic tensor through the two invariant structure 
functions Wi and W2: 



-Wi{prq,Q'){9 



+ 



Wiiprq,Q' 



Pj - Q' 



Pi -q 



u Pj ■ q 

q2 



where j = A,N. Multiplying the left and right sides of Eq. (H) by = — q^ fci_/g_ 
(see e.g. [|T^ ^DQ), where = ei — kiz and g_ = qo — qz, and considering the limit of 
ei, /ci — > 00 with both Q'^ and go fixed, one obtains 



N •' 



2 W) 



— (1+COS5) {p-+q--^ 



\2 I pL ■ 2 X 



(9) 



where sm6 = \/Q^/\q\, P- = Po — Pz, v' = p ■ q/M and w"^ = {p + g)^. Because of the 
off-shellness of the interacting nucleon one has p^ 7^ M^. Contracting Eq. (^) with the unit 
vector = (0, 0, 0, Uy) one has 



N I 



(10) 



''Note that fcf automatically fulfills the current conservation q ■ Jj = and = 0. 
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The inclusive cross section (|]) can be expressed through the structure function and 



Wf in a standard way: 

d^n^ = ^Mou i^W.^iQ', v) + 2tan2(|) W,\Q\ v) 

In case of DIS one introduces the usual scaling functions: 

F{ = MWi 
Fi = uWi, 

where j = A, N. Using Eq. ([I2D in Eqs. (P, p^ , pTD one gets: 



da 



dEf,idVLpi 



cTMott- \f^{x,Q') + ^tan2(|)Fi^(x,g2) 



V 



M 



(12) 



(13) 



where x = tht and 



Fi{x,Q') 



(14) 



— (l + cos5)2(p_ + g_^)2 



pi 



2M2 



sin (5 



(15) 



where i> = (^2 + q2 _ m^)/2M = u' + {p^ - M^)/2M and x = Q^/2Mi>. Note that with 
such a definition of the argument of Fj{x, Q"^) it is ensured that the cross section is vanishing 
below the threshold for the eD e'pn reaction. 

2.1.1 The Nuclear Structure Function F^{x,Q'^) 

In this subsection we will discuss the DIS structure function F^{x, Q"^) in more detail. Let's 
introduce the scaling variables 

Ap_ 



Ma ' 

m7' 



(16) 



and make use of the identity d'^p = ^dp+dp_d'^p±; then, Eq. ([15|) can be written as 
F2^{x, Q') = \Y.j dp^dp^d'p^dz SM F^{i, Q')^5iz - ^) 



N 



■ Ma. 
■AM' 



[1 + cos(5)^(z + a, 



Mv' 



>2 I P± •2 c- 

' ' -sm 



+ 



2M2 



(17) 
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The integration over dp- can be taken automatically, while the integration over which 
describes the virtuality of the interacting nucleon, requires the knowledge of the invariant 
nucleon Spectral Function. One can proceed however by observing that the virtuality of the 
interacting nucleon depends on the structure of the recoil {A — l)-nucleon system. Namely 
for the case of two-body break-up the invariant Spectral Function contains the S{p+ — p+o) 
function with 



p+o = Ma 



{A-z)Ma/A' 



where Mj^_i is the mass of the recoiling A — 1 nucleus. In case of the excitation of the recoil 
nuclear system into its continuum, one can use the observation, based on the multinucleon 
correlation model [^, that for different ranges of z the dominant value of p+ in Sn{p) 
depends on whether the interacting nucleon is in the nuclear mean field or in 2N, 3N, etc. 
correlations. Based on this model we can estimate the integrand in Eq. (|17D as 



]\^2 _j_ 2 

< p+ >= Ma - Ma-2 - -77: Z;-^ at z> 1.2 ^ 1.3 {2N correlations) 

(2 — z)Ma/ A 

(2M)'^ + 

<p+>= Ma- Ma-3 - 7^ — ^ at z > 1.7 ^ 1.8 (3N correlations) 

(3 — z)Ma/A 



(19) 



Using these approximations we can now integrate Eq. ( pTD over p+ arriving at: 

F2^(x,Q2) = ^1 dz£p^ Pn{z,P±) F,^i< X >,g2^ 



where 



< X > 



< u > 



<v > 

'Ma. 21,^ j:.2f ^ M<u' > 2 , pI . 2. 
'AM^ (1 + COS5) (- + «,^^) +^sm5 

_ 

2M<i>>' 
2M 



(20) 



> = + P-1+^ = \P+^i + ^+^] ' (21) 

and p+ here defined according to Eqs. (|18|, [l^). In Eq. (|20|) Pn{z,p±) is the one-body 
density function in the nucleus, defined as: 

Pn{z,p±) = ^ J dp.dp+ Sn{po,Pz,P±)S{z -j^-jj) 
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— j dpo Sn{Po,Po- z—,p±). (22) 



2.1.2 The Bjorken Limit 



Equation ( PU[ ) allows to calculate the inelastic A{e, e')x reaction in a wide range of values of 
Q^, i.e. large enough that the condition for the closure over final hadronic states is achieved 
and the impulse approximation is valid. Additionally, in DIS the range of the Bjorken x 
should correspond to the valence region (x > 0.2^0.3) where shadowing effects are negligible. 
In the Bjorken limit, where Q^, g — > oo and x is kept fixed, Eq. (|l^) transforms to the usual 
convolution formula used by many authors p5|, pOl p2|, E^, p5|]: 



F2^{x,Q') =Y. Tdzz f{z) F^C-.Q% (23) 

A''=l ^ 

where {z) is the nucleon light-cone momentum distribution in the nucleus 

r(z)= f d'p± PNiz,p^), (24) 



with the baryon charge normalization condition given by Eq. (pTj). 

Introducing the compact notation f'^ to indicate the convolution (|23D and assuming 



exact nuclear charge symmetry, the ^He and DIS structure functions can be written as 

^2^^= = S^{2Fi + F^) + D^{2FP-F^), 

= S®i2F^ + FP) + D®{2F^-FP), (25) 



where 



S{z) 



p{z) + r{z) 



2 



with {z) being the light-cone momentum distribution of proton (neutron) in ^He. If 
D{z) ~ [i.e. P{z) ~ /"(-z)], then it is reasonable to expect that the EMC ratios (|I|) 
in ^He and are quite close each other, so that the super-ratio (g) is close to unity, as 
observed in Ref. p5[. However, as it will be illustrated in detail in the next subsection, the 



spin-flavor dependence of the nuclear force (even without any charge-symmetry and charge- 
independence breaking terms) yield f^{z) ^ f"'{z). Therefore, when D{z) ^ 0, the difference 
in the proton and neutron structure function [leading to 2F2 — F2 7^ 2F2" — in Eq. (P3|)] 
can give rise to RJmc Remc correspondingly to deviations of the super-ratio (H) from 
unity depending on the size of the EMC effect itself. It is important to note that the nuclear 
charge symmetry will not limit such deviations. 
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2.2 Nuclear Density Function and LC Momentum Distribution in 
the VNC Model 

Now we further analyze the one-body density function Pn{z,P±) and the hght-cone momen- 
tum distribution {z). From Eqs. (0) and (^) one obtains the following normalization for 
the one-body density function [pO[] : 



Jo 



A pA 

dz z PN{z-iPi_) = dz z f^{z) = 1. (27) 
Jo 



To construct the one-body density function and subsequently the light-cone momentum 
distribution we have to relate the invariant Spectral Function Sn{p) to the non-relativistic 
Spectral Function P'^{p,E), which represents the joint probability to find in the nucleus a 
nucleon with three-momentum p = \p\ and removal energy E. Since the latter is defined as 
E = Ea—Ea-i+E\_-^ [Ea-1 being the (positive) excitation energy of the system with (A — l) 
nucleons measured with respect to its ground-state, and Ea {Ea~i) the binding energy of 
the nucleus A {A — 1)], the nucleon Spectral Function also represents the probability that, 
after a nucleon with momentum p is removed from the target, the residual {A — l)-nucleon 
system is left with excitation energy E\_^. 

Since the (non-relativistic) nucleon Spectral Function P^{p,E) is normalized as 



Att dE dp p' {p, E) = 1, (28) 

where Emin = Eia—Ea-\ is the minimum value of the removal energy, one has some ambiguity 
in the relation between Sn{p) and P^{p,E). Two ansatz were suggested to relate Sn{p) 
and P^ip^E) which can be considered to represent two extremes. In one [17, 30, 31] it is 
assumed that 

Sn{p) = ^-P\p,E). (29) 
Apo 

In this case the non-relativistic transition from Sn{p) to P^{p,E) is straightforward, since 
in this limit one has {Ma/Apq) — >-~ 1, and hence the renormalization is the smallest for 
small nucleon momenta. Another ansatz assumes that renormalization is momentum 
independent so that 

Sn{p)^Cn-Pn{p,E), (30) 

where can be found from the requirement given by Eq. (|^). In this prescription f^{z) 
reads explicitly as 

f{z) = 27,MCM dE dppP''{p,E)^==^, (31) 



where Pmin{z, E) is given by 



p„.„(..E)a(^-f)'-'-'^-.)'. (32) 

i"mtn\ '7 2 Ma-Mz ^ ' 
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with Ma = a ■ M + Ea and M^_i = Ma~i + E — Emin- Note that the normahzation factor 
Cm can be different for protons and neutrons, in hne with the normahzation factor appearing 
in Eq. 

For a generic nucleus with A > 2 the nucleon Spectral Function can be written as (cf., 



e.g., H, 

P^(p, E) = Po^(p, E) + (p, E), (33) 

where P^{p, E) includes the contributions of all the final states belonging to the discrete 
spectrum of the {A — l)-nucleon system (basically its ground and one- hole states), while 
P^{p,E) corresponds to more complex final configurations (i.e. the final states of the con- 
tinuum spectrum of the [A — l)-nucleon system), which are mainly lp-2h states arising from 
the 2p-2h excitations generated in the target ground state by short-range and tensor NN 
correlations. In what follows we will refer to and P{^ as the ground and correlated parts 
of the nucleon Spectral Function, respectively. 

The nucleon momentum distribution n^{p) can be simply obtained from the nucleon 
Spectral Function by integrating over the removal energy; thus, Eq. (^) implies that n^{p) 
can be written as the sum of two components related to the ground and correlated parts of 
P^{p,E), respectively, viz. 



(p) ^ / dE P^(p, E) = n^ip) + nUp). (34) 



Useful parameterizations of the results of many-body calculations of n'^{p), available for 
few-nucleon systems, complex nuclei and nuclear matter, as well as its decomposition ( p4D 
into Uq (p) and n^{p) can be read off from Ref. P3|(b). 

As it is well known, the calculation of P'^lp, E) for A > 2 requires the knowledge of 
a complete set of wave functions for {A — 1) interacting nucleons. Thus, since the latter 
ones should be obtained from many-body calculations using realistic models of the NN 
interaction, the evaluation of P^{j),E) represents a formidable task. In case of ^He the 



nucleon Spectral Function has been obtained using three-body Faddeev |36] or variational 
37| wave functions, whereas for A = oo the evaluation of P^ {p, E) has been performed using 



the orthogonal correlated basis approach |Q and perturbation expansions of the one-nucleon 



propagator Since in this work we are interested in the evaluation of Eqs. (p0|j23|) for 
mirror A = 3 nuclei as well as for A > ?> nuclei, in what follows we will adopt the Spectral 
Function model of Ref. , which was developed for any value of A and shown to reproduce 
in a very satisfactory way the nucleon Spectral Function in ^He and nuclear matter calculated 
within many-body approaches using realistic models of the A^A^ interaction. 

Let us now briefly describe the (non-relativistic) nucleon Spectral Function adopted in 
case of '^He and ^H. For the former nucleus the ground component P^{p, E) is given by 



mm 



PSip,E) = nl{p)5[E-E, 
P^{p,E) = 0, (35) 
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where UqIp) is the proton momentum distribution corresponding to the ^He to deuteron 

transition and eU^^"* ~ 5.49 MeV . Note that the ground component Pq{p-, E) is identically 
vanishing because the residual pp system does not possess any bound states. In case of ^if, 
since charge symmetry largely holds for the nuclear wave functions of ^He and ^if, one can 
write: 

P^{p,E) = R{p) . nl{p) 5[E - EifJ] , 

P^{p,E) = 0, (36) 



where n^{p) is the same momentum distribution appearing in Eq. (pSl) . The correction 
function -R(p), which includes the effects due to charge symmetry and charge independence 
breaking terms in the A^A^ (as well as NNN) interaction, turns out to be quite close to 
unity, namely within 2 -h 3% level of accuracy, as it can be estimated from the explicit 
calculations of ^He and wave functions carried out in Ref. For n^{p) we use the 



simple parameterization obtained in Ref. [^](b) in case of the RSC model |0 of the NN 

interaction. Note that -E^in — 6.26 MeV ^ E)^^^ because of the different values of the 
experimental binding energies of ^He and . 

As for the correlated part Pf{p,E) for ^He we adopt the model of Ref. |3^ and the 



parameterizations of "'^(p) corresponding to the RSC interaction. For we follow 



the same logic of Eq. (|36|) and estimate the correction function R{p) from the calculations 
of Ref. [0. An additional difference in the correlated parts may arise from the different 
values of the threshold for the three-body break-up Ethr (cf. [§3|)' namely: Ethr = 7.72 and 
8.48 MeV in ^He and ^H, respectively. 

In Fig. 1 we have reported the results for two ratios of light-cone momentum distribution 
functions, namely f^He/f^H (solid line) and f^Hel f^u (dashed line), estimated according to 
the results of Ref. |^ . It can clearly be seen that the corrections to the charge independence 
in case of ^He and nuclear wave functions are typically at the level of 2 3%. 

Next we want to estimate the uncertainty introduced by the above-mentioned normaliza- 
tion procedures, given by Eqs. (^) and (|30D . In Fig. 2 we have reported the results of our 
calculation of the proton and neutron light-cone momentum distributions in ^He according 
to Eqs. (|30D and (|3TD (solid and dashed curves for proton and neutron, respectively) and 
according to Eqs. (^4]) and ( ^Of ) (solid and dashed curves with circles for proton and neutron, 
respectively). It turns out that the two different normalization prescriptions can substan- 
tially differ at small values of the light-cone fraction z and this may represent a potential 
source of uncertainty. 

From Fig. 2 it can be seen also that within each normalization prescription /^(z) 7^ f^{z). 
Such a difference is driven by the presence in the three-nucleon wave function of a mixed- 
symmetry S"-wave component as well as of P- and D-waves arising from the spin-spin, 
spin-orbit and tensor terms of the A^A^ interaction, respectively. According to the results 
of sophisticated solutions of the three-nucleon ground states both with and without charge 



symmetry (and charge independence) breaking terms |4l|, the probabilities of the S", P 



and D partial waves are: Ps' — 1.2 ^ 1.5%, Pp < 0.2% and Pd — 7 ^ 9%, depending on 
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the specific model adopted for tlie NN (as well as NNN) interaction. Therefore, we do 
not expect that the dependence of f^{z) on the particular nuclear force model could be 
significant for the estimate of the deviation of the super-ratio @ from unity, as it is also 



suggested by the results already obtained in Ref. [g^. Finally, we want to stress that our 
results for the EMC ratio (|l]) obtained within the convolution formula (p3D for mirror A = 3 
nuclei, have been positively checked against the corresponding results of Ref. and [p6| , 
where the nucleon Spectral Function obtained from few-body variational techniques in case 
of the RSC potential has been employed; we have found that the differences do not exceed 
~ 1% in the whole x-range of interest in this work, viz. 0.3 < a; < 0.9. 



2.3 Light-Cone Approach 

In this subsection we discuss the formalism of LC dynamics in the approximation where the 
relevant degrees of freedom in the nuclear medium are the nucleons only, carrying therefore 
the total momentum of the nucleus. The structure function in the LC framework can be 



represented as follows |T^, 



Af=l ^ 



< U > 



^^'{1 + cos5)2(z + a /" ^'^ 



'AM' 



M < z/' > , o . 9 ' 



2M2 



(37) 
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Figure 2. The nucleon light-cone momentum distribution in ^He, corresponding to the RSC 
model |4C] of the NN interaction, adopted in this work. The dashed and solid lines correspond to 
the neutron and proton momentum distribution, respectively. Lines with open circles correspond 
to the calculation with the normalization scheme of Eq. (p9|) , and those without open circles to the 
normalization scheme of Eq. 



where all the quantities in the r.h.s, except p'ff^ZjP^), are defined in subsection In 
the Bjorken limit one obtains: 

F^^i^, Q') = E r -/^"^ W i^2^(-, Q% (38) 

N=l •'^ ^ ^ 

where fjf^{z) is related to Pff{z,p±) according to Eq. (|2^) . The quantity pjf{z,p±) rep- 
resents the nucleon LC density matrix in the nuclear medium. This function satisfies two 
sum rules; namely, from baryon charge conservation one has 

|^dVp^^(^,P±) = l, (39) 
while the momentum sum rule requires 

[-d'p^zp^^^{z,p^) = l. (40) 

Note that the last sum rule is not directly satisfied in the VNC model, but it can be restored 
if mesonic degrees of freedom are explicitly introduced. 

In general p^{z,pA_) is not known for nuclei with A > 3. However, for numerical calcu- 
lations one can proceed using the two following observations. First, in the non-relativistic 
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limit (applicable for 0.7 ^ 0.8 < ^ < 1.2 ^ 1.3 and p±<kF^ 200 ^ 300MeV/c) the density 
p^{z,p±) can be related to the non-relativistic nucleon momentum distribution n^{p) as 



N/ 



(41) 



Second, within the two-nucleon correlation model one can relate the high momentum 
tail of the nuclear LC density matrix to the two-nucleon density matrix on the light cone 



0: 



Pn''{z,P±) 



02(^4) p^^^{z,pa_) 



Ek 02 (A) riNNik) 



(42) 



where Ek = VM^ + and p 



m'^+pI 



-M2. 



(43) 



The observation that the two-nucleon correlations define the high momentum tail of the 
nuclear momentum distribution at momenta (> kp) allows to replace n^{k) in Eq. (|4^) by 
a2{A) nNN{k) at > kp, where a2{A) characterizes the probability to find a two-nucleon 
correlation in the high momentum tail of the nucleon momentum distribution in the nucleus 
y40. Finally in the kinematical range where two-nucleon correlations dominate the proton- 
neutron pair, one can replace riNNik) by the squared LC wave function |^'£)(/c)p of the 
deuteron P^ . 

The LC many-nucleon approximation for the nuclear wave function leads to a prediction 
for the F2IF2 ratio which qualitatively contradicts the EMC effect for x > 0.5. This 
reflects the need to include explicitly non-nucleonic degrees of freedom in nuclei in order to 
explain the EMC effect. In the LC approximation a natural explanation is offered by the 
deformation of the quark wave function in the bound nucleon which will be considered in the 



subsection [4.4| . In the following Section we will use the LC model to illustrate the magnitude 
of the Fermi motion effects on the super-ratio (0). 



3 Numerical Estimates of Inelastic Cross Section 



To check the reliability of our assumptions in the derivation of Eq. (pOD as well as of the 
models used for the three-nucleon Spectral Function, we first compare our calculations with 
the experimental data at moderate values of ~ 2 ^3 {GeV jc)^ . At these values of the 
nuclear modification of the valence quark distributions (the EMC effect) is expected to be 
small and thus the comparison with the data will allow to check the validity of Eqs. (pO|-p^). 



Figure 3 presents the comparison of our calculations with the experimental data of Ref. [^, 
where for the nucleon structure function F^ix^Q"^^ we have used the parameterization of 

^An estimate of 02 (A) for a variety of nuclei can be found in Ref. Q. 
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Ref. containing also the contribution of nucleon resonances. For the evaluation of [see 
Eq. dlD] we have used the relation = Fi^[l + (2Mx/z/)]/2x(l + R) with R = 0.18. The 
comparison clearly demonstrates that Eq. (^) is a good starting point for the discussion of 
higher regime. 
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Figure 3. The cross section of inclusive ^He{e, e')X scattering as a function of energy transfer u. 
Dashed hne - inelastic contribution calculated according to Eq. ( pO| ) and adopting for the nucleon 
structure function the parameterization of Ref. |^], which includes nucleon resonances. Dotted line 
- quasi-elastic contribution calculated according to Ref. |44|. Solid line - total cross section. The 
experimental data (full triangles) are from Ref. [^ ] and the kinematical conditions are shown in 
the insets. 

At larger the first question we want to address is how fast the Bjorken limit is 
established and how much the nuclear recoil effects accounted for in Eq. (^) are important. 
The finite effects are governed by the scale of the target-mass corrections (~ M'^/Q'^ as 
well as the factors proportional to Q^jv^ and p\/Q'^)- In Fig. 4 we compare the results 
obtained for both the EMC ratio and the super-ratio for '^He and targets at = 
10 (GeV/c)'^, calculated within Eq. (pO|) and the convolution formula ([23|). It can clearly be 
seen that, while for the EMC ratio (P the convolution formula works within a 2 3% level 
of accuracy at large x, in case of the super-ratio (^ the differences between the non-Bjorken 
and the Bjorken limits cancel out almost completely. 

The next question is the expected uncertainty on the EMC ratio due to the different 
normalization procedures of the nucleon Spectral Function discussed in the previous Section. 
In Fig. 5 we compare the calculation of TIeIac ^^"^ STZemc performed within the convolution 
approximation using the two different schemes of normalization given by Eqs. 



^ and (|3g) 

It can be seen that "R-e^aw exhibits some sensitivity to the choice of the normalization scheme 
while the differences in the calculated STZemc are well below ~ 1%. 
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Figure 4. The x dependence of T^emc (^) STZemc (b) for ^He and targets at = 
10 {GeV/c)"^. Dashed and sohd hnes are the results obtained using Eq. (|20|) calculated without 



invoking the Bjorken limit, and Eq. (23) using the Bjorken limit, respectively. Lines marked by 
crosses correspond to target, unmarked lines to ^He target. In (a) and (b) the CTEQ set 



of PDF^s from Ref. [10| is adopted, while in (b) the results obtained using the parameterization 
of Ref. for the nucleon structure function F^{x,Q'^) are also reported. The charge-symmetry 
breaking effects shown in Fig. 1 are included in the calculations. 



Next we address the sensitivity of the super-ratio (Q) to the particular choice of the 
PDF parameterization in the nucleon. To this end we have calculated the super-ratio (§) 
via the convolution formula (|23| ) using the proton and neutron LC momentum distributions 



(pT|) with the normalization procedure given by Eq. (pOl). The RSC model of the A^A^ 
interaction adopted in the calculation yields Cn = 1.048 and Cp = 1.033 for the normalization 
constants in Eq. (|30|) . We have neglected the charge-symmetry breaking effects shown in Fig. 
1 and we have used different parameterizations of the nucleon structure function F^{x, Q^) 
taken at = iq {GeV/cf, namely the GRV set Q of PDF's and the SLAC fit of Ref. 
1^. Both the GRV and SLAC parameterizations are constructed in a such a way that the 
neutron to proton structure function ratio reaches the "non-perturbative prediction" 1/4 as 
X 1 (see Fig. 6(a)). We have therefore applied to the GRV and SLAC structure functions 
an ad hoc modification in the form of a distortion of the c?-quark distribution limited only at 
large x (i.e., x > 0.7), viz.: d{x) d{x) + 0. 1x^^(1 + x)u{x). Such a modification has been 
directly implemented in the GRV set of PDF^s, while in case of the SLAC parameterization 
we have considered the following replacements: -Ff (x, Q^) Pii^y + 0. 1x^^(1 + x)/4} 

and F2{x, Q^) — >■ F2{x, Q^){1 + 4 ■ 0.1x'^(l + x)}. In both cases the n/p ratio of the modified 
structure functions goes to the "pQCD prediction" 3/7 ~ 0.43 as x — 1 (see Fig. 6(a)). 
Since the proton structure function is dominated at large x by the w-quark distribution, the 
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Figure 5. (a) The EMC ratio [Eq. ^] in ^He (thin hnes) and in (thick hnes) vs. x at 
= 10 (GeV/c)'^. (b) The super-ratio [Eq. (|2|)] vs. x. Dashed and solid lines correspond to the 
results obtained using the normalization schemes of Eqs. (^) and (|30|), respectively. The CTEQ 
set of PDF^s from Ref. |1C] has been adopted as input for the nucleon structure function F2 {x, Q^). 
The charge-symmetry breaking effects shown in Fig. 1 are not included in the calculations. 



above-mentioned modification does not change significantly the shape of the proton structure 
function; the effects are larger on the neutron structure function, but by construction they 
are limited in the region x > 0.7 (see Fig. 6(a)). 

The results obtained for the super-ratio @ are shown in Fig. 6(b). It can be clearly seen 
that the deviation of the super-ratio from unity is small (less than 1%) up to a; ~ 0.75, while 
it increases rapidly as a; > 0.75 and depends strongly on the large-x behavior of the n/p 
ratio. Our conclusion is that the VNC model predicts a deviation of the super-ratio from 



unity within 1% only for x < 0.75 in overall agreement with the results of Refs. p5|, ^ 
Note however that the x-shape and the average value of our results for the super-ratio are 
closer to the findings of Ref. |26| (where a Spectral Function similar to the one of the present 
work is adopted) and differs from the results of Ref. where larger deviations (up to 

2%) from unity were found. It is likely that the difference is related to the different Spectral 



Functions used in the present work and in Ref. ||25|, since the latter uses the VNC model 



with the same normalization scheme [Eq. ([29|))] for the nucleon LC momentum distribution. 

The x-shape of the super-ratio S71emc{.x,Q'^) shown in Fig. 6(b), can be better under- 
stood by looking at Fig. 7, where the EMC ratio 'R-emc{^iQ'^) iii the two mirror A = ?> 
nuclei is separately reported. It can be seen that the convolution approach predicts a larger 
deviation from unity in than in ^He for x < 0.75. This is a direct consequence of the 
higher kinetic energy of the neutron (proton) with respect to the proton (neutron) in ^He 
{^H) due to the spin- flavor dependence of the nuclear force [see the discussion after Eq. (P5|)]. 



18 



o 

>< 



0.8 



0.6 



0.4 - 



o 

'Lj.'^ 0.2 



T — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 




(a) 



Q Q I I I I I I I I I I I I I I I I I I I I I I I I I I I u 



O 



1 .02 



1 .00 



0.98 - 



0.96 - 



0.94 



0.92 - 



0.90 



T — I — I — I — I — I — I — I — I — I — I — |-T — I — I — I — I — I — I — I — I — I — I — I — I — I — r 



(b) 



J I I I I I I I I I I l_l I L 




0.3 



0.5 



0.7 



0.9 



0.3 



0.5 



0.7 



0.9 



X X 

Figure 6. (a) Neutron to proton structure function ratio, ^2^(2;, Q^)/F2 (x, Q^), vs. the Bjorken 
variable x at = 10 (GeV/c)'^. The full dots are the NMC data points as given in Ref. p5|. Thick 



dashed and solid lines correspond to the GRV set |11] of PDF^s and to the SLAC parameterization 
of Ref. [^. Thin lines are the modified GRV and SLAC fits, as described in the text, (b) Super- 
ratio [Eq. (|2|)] of the EMC effects in yl = 3 mirror nuclei. The meaning of the lines is the same as 



in (a). Using the CTEQ parameterization |10| one obtains results very similar to those reported 
for the GRV set llTTI of RDTs. 



The super-ratio S71emc{x,Q'^) is therefore larger than one up to x ~ 0.75 (see Fig. 6(b)), 
but such a deviation from unity is small because the EMC ratio itself is predicted to be 
quite small in the two mirror nuclei (less than a 1% effect) within the VNCM. For x > 0.75 
the EMC ratio increases above unity very sharply; generally speaking, this is related to 
the fact that the nucleon structure function goes to zero as a; — > 1, while the nuclear one 
is non-vanishing because of the Fermi motion of the nucleons in the nucleus. Moreover, the 
slope of the rise of TZemc{x, Q^) is larger in than in ^He due to the decrease of F^/Ff 
at X — > 1. Note also that the EMC ratio in is sensitive to the modification of the d/u 
ratio at large x, whereas the EMC ratio in ^He is not (see Fig. 7). Thus, for x > 0.75 
the super-ratio STlEMci^yQ^) drops below one and becomes a rapidly varying function of 
X with a remarkable sensitivity to the large-x shape of the nucleon structure function. 

Finally the last source of uncertainty we want to consider within the convolution approx- 
imations is the difference between the predictions of the VNC model and the LC formalism. 
In Fig. 8 the predictions for EMC ratio calculated within the VNC model according to Eq. 
(pop and the normalization scheme of Eq. (|29| ) are compared with the corresponding ones 
of the LC approach [see Eq. (|37D]. In both cases we have adopted the parameterization 
of Ref. [^] , which contains the contribution of nucleon resonances. As it follows from Fig. 
8(a) the LC approximation predicts larger value of the EMC ratio as compared with the 
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Figure 7. The EMC ratio [Eq. (|l])] in ^He (thin lines) and in (thick hnes) as predicted by 
the convolution formula ( |2^ ) at = 10 (GeV/c)'^. The solid lines corr espond to the SLAC || 
parameterization of the nucleon structure function, while the dashed lines are the results obtained 
using the modified SLAC fit at large x as described in the text. The thin dashed and solid lines 
are almost indistinguishable. 

VNC model. As a result the super-ratio within the LC approximation is smaller (closer 
to one) as compared with the prediction of the VNC model. Note also that the effects of 
nucleon resonances are still visible in Fig. 8(a) for x > 0.8, corresponding to < 2 GeV at 
= 10 {GeV/cy. Therefore, if one wants to investigate only the leading twist of the nu- 
cleon structure function, one can either limit the range of values of x or increase sufficiently 
the value of Q^. 

To sum up this Section, we conclude that all the considered uncertainties within the 
convolution approximation, in which no nuclear modification of bound nucleons is considered, 
do not yield deviations of the super-ratio (^ from unity larger than 1% at a; < 0.75 (2% at 
x<0.8). 

4 Models of the EMC Effect with Modifications of the 
Bound Nucleon Structure Function 

Although within the VNC model and the LC approach the nuclear corrections to the super- 
ratio (^ are ~ 1% at x < 0.75, it is hardly safe to treat this as an ultimate estimate of the 
nuclear effects. The VNC model is just one of the many models of the EMC effect. Also, 
literally, the VNC model predicts the parton densities to violate the momentum sum rule 
(for instance by ~ 5% for an iron target). When this feature is fixed by adding mesonic 
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Figure 8. (a) The EMC ratio [Eq. ^] in ^He (thin hnes) and in (thick hnes) vs. x at 
= 10 (GeV/c)'^. (b) The super-ratio [Eq. (|2|)] vs. x. Dashed and sohd lines correspond to the 
results obtained using the VNC model and the LC formalism, respectively. The parameterization 
of Ref. which includes nucleon resonances, has been adopted as input for the nucleon structure 
function F2 [x, Q"^)- The charge-symmetry breaking effects shown in Fig. 1 are not included in the 
calculations. 



(pionic) degrees of freedom, one predicts an enhancement of the antiquark distributions in 
nuclei at x > 0.05 which grossly contradicts the Drell-Yan data [Q. It is also well known 
that the convolution approximations underestimate significantly the EMC effect at large x 
(cf., e.g., Refs. |T^, ^ and |^). Experimental data are available for a variety of nuclei and 



in Fig. 9 we have limited ourselves to the cases of ^He and ^^Fe. The convolution formula 
within the VNC model ( P3| ) has been evaluated adopting for the nucleon Spectral Function 
P'^lk, E) the model of Ref. and our results turn out to be in agreement with the findings 
of Ref. From Fig. 9 it can clearly be seen that the convolution approach is not able to 
reproduce the minimum of the EMC ratio around x ~ 0.7 as well as the subsequent sharp 
rise at larger x. Note that the disagreement is even larger within the LC approximation (see 
dashed curves in Fig. 16). 

Therefore, it is reasonable to expect that the results of the convolution approximation 
for mirror A = 3 nuclei can suffer the same drawback. Moreover it is very important to 
asses any isospin dependence of EMC effect in order to extract in a reliable way the neutron 
structure function from ^He and data. An isospin dependence for EMC effect is naturally 
expected from the differences in the relative motion of pn and nn {pp) pairs in i^He). 
The results obtained for in Ref. in case of the Argonne V"18 + Urbana IX models 
of the NN and NNN interactions, are reported in Fig. 10. It can clearly be seen that, 
since the interaction of a pn pair is more attractive than the one of a nn pair, the proton 
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Figure 9. The EMC ratio [Eq. (|)] in ^He (a) and ^^Fe (b) at = iq [QeV/cf. Open dots 
are data from Ref. [46|, while in (b) the full squares are from Ref. ||4^. The solid lines are the 
results of the convolution formula (|2^), calculated adopting the SLAC |^] parameterization of the 
nucleon structure function F2 {x, Q^) and the model of Ref. [33| for the nucleon Spectral Function 
P^{k,E). 



is closer to the center-of-mass than the neutron. The corresponding root mean square 



radius turns out to be: ^< r^^ > ~ 2.5 fm and y< r^^ > ^ 2.8 fm. As a consequence, 
the overlapping probability is larger for a pn pair than for a nn pair. As a matter of fact, 
from Fig. 10 it follows that the partially integrated probability to find a A^A^ pair with 
^NN ^ ^ f^, is ~ 40% larger for a pn pair than for a nn pair. We stress that this is a very 
important isospin effect in mirror A = 3 nuclei. 

Thus in order to draw final conclusion about the size of the deviation of the super-ratio 
<STZemc{x,Q'^) from unity we should investigate effects beyond those predicted by the convo- 
lution approach. 

The discovery of the EMC effect at large x has triggered a huge theoretical effort which 



has led to the development of a large number of models (see, e.g., Refs. |^9|, |T6|, ^ 
and references therein). In this work we will limit ourselves to consider some of these 
models, which are of interest for an estimate of the possible deviation of the super-ratio 
STZEMci^y Q^) from unity. We will use the experimental points presented in Fig. 9 in order 
to constrain as much as possible the parameters of these models. Note that the experimental 
uncertainties on the EMC ratio in ^ife are significantly larger than the corresponding ones 
in ^^Fe; therefore, new measurements on '^He targets with reduced errors will certainly help 
in improving our knowledge of the EMC effect in light nuclei. 
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Figure 10. The distributions gNNii^NN) of the relative motion of a nn pair (fuU dots) and of a 
pn pair (open squares) in ^H, as a function of the relative distance ttvat between the A^A^ pair. 
The results reported correspond to the the Argonne V18 and Urbana IX models of the A^A^ and 
NNN interactions, respectively, obtained using the Green Function Monte Carlo method of Ref. 
[42|. The distributions are normalized as: drj\f]\fgj\f]\f{r]\fN) = 1. 



4.1 Nuclear Density Model 

It was argued in [T^, ^ that due to a diluteness of the nuclear systems the nuclear effects 
for the deviation of the nuclear structure function from the sum of the nucleon structure 
functions can be treated as a series in the powers of /c^/M^ and e^/M. This approximation 
holds in a number of dynamical models, like the rescaling model f^ , the six quark 

(6g) cluster model ||51|, |52|, the color screening model (suppression of small size configurations 



in bound nucleons) ||T^, ^ and pion models |Q . Hence in the region of small enough x 
(i.e, X < 0.7), where terms oc k'^/M'^ can be neglected, an approximate factorization should 
take place 



T^EMc{^. Q') - 1 = Q'') fiA), (44) 
where TZiMci^^ Q^) = F^i^, Q^)/[ZF^{x, Q^) + NF^{x, Q^)], and 

f{A) oc< k^ > /M^ (45) 



or to the average virtuality of the nucleon. Eq. ( ^^ is in a very good agreement with 
the SLAC data on the A-dependence of the EMC effect. Numerical estimates using Eqs. 
(K3,^) and realistic deuteron and iron wave functions lead to 



Ff (x, Q2) + F^{x, Q2) 4 (x, Q^) 
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For A > 12 one has approximately < k"^ > /M^ oc< Pa(^) >= p(^)) where p{A) is the 
average nuclear matter density, leading to 



Kmc{^, Q') ^ ^r^j^^'l; = a{x, Q') [1 + p{A)P{x, Q')] . (47) 



Note in passing that such an approximation is definitely not applicable at very large x, since 
short-range correlations dominate for a; > 1 and therefore the relation TZ^^j^^x) oc pA is 
expected not to hold any more. Also one hardly can directly use this approximation for the 
deuteron since the notion of average nuclear density is not well defined in this case. 

Analysis of the data on the EMC effect using Eq. (^) (including deuteron and "^He 
data) has been carried out in terms of the average nuclear density p{A) in Ref. [^]. The 
quantities a(x, Q^) and P{x,Q^) were fitted to the data; their values for various x-bins can 
be easily read off from Table IX of [HH]. The nuclear density p{A) was assumed to be given 



by : p{A) = 3A/4:7iRl where Re = yS/S ■ r^^A, with r^^^ representing the r.m.s. electron 
scattering (charge) radius of the nucleus. In particular, the values pi^He) = 0.089 fm~^ 
and pi^^Fe) = 0.117 fm^^ were adopted in We will refer hereafter to Eq. ( ^7|) as the 



density model. 

Note that the fit in Eq. (^T]) has been done in Ref. using charge rather than matter 



radii of nuclei, which is a good approximation for large A and Z = N nuclei since in this 
case 

<^ ^ = <^ ^ -I- ^ ^ -I- ^ ^ { 48l 

c,A ^ matter, A ^ ~ ^ ' ^proton ^ ~ ^ ' c,neutron \ "/ 



and 



'^matter, A ^ ^ ^c,proton ^ ' ^ ^c,neutron ^ ■ (^'^) 



For light isosinglet nuclei Eq. (^) is expected to hold. However the predictions of the 
density model for light isosinglet nuclei should have a rather qualitative character since Eq. 
(^) does not hold for deuteron and barely holds for ^He nucleus. Moreover, the step leading 
from Eqs. (||J4|) to Eq. (|^ is not justified. 

We have mentioned above that the density model was proposed in Ref. [|l6l only in 
case of sufficiently heavy nuclei. If we want to apply Eq. ( ^Tj ) to mirror A = 3 nuclei, the 
first question is which density we have to use. As already observed in Fig. 10, the neutron 
(proton) is closer to the ^He {^H) center-of-mass than the proton (neutron). This means 
that the neutron (proton) has more kinetic energy of the proton (neutron) in ^He {^H)^ 
According to the RSC interaction, the neutron in ^He possesses on average about 25% 
kinetic energy more than the proton. Since the deviation of the EMC ratio from unity may 
be related to the mean kinetic energy of the nucleon and to the derivatives of the nucleon 
structure function (cf. |T^), we expect a different EMC effect in mirror A = 3 nuclei, driven 

''in what follows we neglect the small isospin violation driven by charge symmetry breaking effects in the 
relative motion of the pp pair in "^iJe and of the nn pair in '^H. 
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by the spin-flavor dependence of the A^A^ interaction and by the different quark content of 
the proton and neutron (cf. Fig. 6). In case of A = 3 systems it should be emphasized that 
what matters ultimately is the matter size and not the charge radius. The relation between 
charge and matter radii for ^ife and targets differs from Eq. (0), namely: 

<^ -i ^ = ^ ^ -I- ^ ^ A- <^ ^ /2 

c,^He matter,proton ^ ~r ^ ' c,proton ^ ~r ^ ' c,neutron ^ / 

^ ^c,^H ^ ^ ^matter, neutron ^ ~l~ ^ ^ c,proton ^ ~l~ 2 < c,neutron ^i (^^) 

where < rl^atter, proton > and < rl^atter,neutron > are matter density of proton and neutron in 
^He. Thus a relation similar to Eq. ( ^8| ) can be obtained only for the isosinglet combination 

of the ^He and targets since in this case < rl^3He+^H)/2 > = < rl^atter,{proton+neutron)/2 > 

+ < "^Xproton > + 5 < rl neutron > /4, which coiucidcs with Eq. (H) up to the small term 

^ ^c,neutron ^ 

In case of ^He and nuclei we need to account for the fact that it is the difference 
between matter radii of proton and neutron in ^He (or in ^H) that should be considered in 
estimating the different EMC effects for these nuclei within the density model. To be able to 
use the results of the fit of Ref . one should use for proton (neutron) density Pp(„) (A) = 

3A/47ri?3(„), where Rp^n) = ■ Vp^n) and r^^n) =< riatter,protonineutron) > + ^Q, wherc the 



parameter tq accounts for the fact that the fit of Ref. is based on the use of the nu- 
clear charge radius. For estimation purposes we take tq ~ ^ < r^^proton > + < fc,neutron > — 
0.7 fm. 

Using these densities one can now estimate the EMC effects within the density model, 
modifying the Eq. ( P7| ) as follows: 

^EMC ~ a i + p 



2F|(x,Q2) + Fn(a;,g2^ 



In Fig. 11 the predictions of the density model ( ^71) for both the EMC ratio and the 
super-ratio are reported and compared with the results of the VNC model. For proton 



and neutron matter radii we adopt the values J< r^atter,proton > = 1-75 ± 0.03 fm and 



< f^rnatter,neutron > = 1-55 ± 0.04 fm, obtained from the results of Refs. ||55|, ^ p7 |. 
Since in this simple model we neglect small effects of the isospin violation which could 
lead to < rl^atter,proton ^ i^He) rl^atter, neutron > i^H), the above rcsults correspond to: 
Ppi^He) = Pni^H) ~ 0.050 /m'^ and Pn{^He) = pp{^H) ~ 0.068 /m'^. 

It can be seen from Fig. 11 that the deviation of the EMC ratios from unity is different 
for ^He and targets by approximately the same amount in percentage, but the density 
model predicts a deeper EMC effect. Therefore, at variance with the VNC model, the 
deviation of the super-ratio STZemc{.x, Q'^) from unity can reach a ~ 2% level already around 
X ~ 0.7 -h 0.8 in the density model, because the latter predicts a larger EMC effect with 
respect to the VNC formula. 
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Figure 11. (a) The EMC ratio [Eq. (||)] in ^He (thin hnes) and in (thick Unes) vs. x at 
= 10 (GeV/c)'^. (b) The super-ratio [Eq. (^)] vs. x. The meaning of the Unes is the same as in 
Fig. 7. In (a) the open and full dots correspond to the predictions of the density model (^) for 
^He and ^H, respectively, adopting ppi^He) = Pni^H) = 0.050 /m~^ and pn(^He) = Ppi^H) = 
0.068 fm^"^. In (b) the prediction of the density model is represented by the open squares. 



However we stress again that one should be very careful in applying the density model for 
light nuclei, as ^He and ^H, since for the lightest nuclei the Fermi momentum distribution 
is very steep and the A-dependence of the EMC effect may not have the same form as the 
one for heavy nuclei. Therefore the predictions of the EMC effect for ^He and targets 
based on the density model should be considered for illustrative purposes only. 



4.2 Quark Confinement Size 



In Refs. it was proposed to explain the EMC effect at large x via the softening of 

the (valence) quark distributions in nuclei (i.e., a more efficient gluon radiation in bound 
nucleons than in free nucleons) caused by an increase of the confinement volume of the quark 
in a bound nucleon. One can combine the model of with the VNC model by including 
in the latter modifications of the structure functions of the virtual nucleons. For simplicity 
one can neglect the dependence of the modification on the nucleon momentum treating this 
effect on average. In this case one can write 



x 



(52) 
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where ^a(Q^) is the (dynamical) rescahng factor, whose Q^-dependence, dictated by pQCD, 
is given by 



X2 



(53) 



with Xa and representing the quark confinement sizes in the bound and free nucleon, 
respectively. A change of Xa with respect to A at may be viewed as a change in the nucleon 
size in the nuclear medium (this interpretation is usually refereed to as the nucleon swelling). 
In this respect, it should be pointed out that: i) in Ref. |Q an increase not larger than 
~ 6% of the proton charge radius is found to be compatible with y-scaling in ^He and ^^Fe; 
ii) the analysis of the Coulomb Sum Rule (CSR) made in Ref. suggests an upper limit of 
~ 10% for the change of the proton charge radius in ^^Fe; iii) recently the experimental 
values of the CSR in ^^C and ^^Fe have been re-analyzed at ~ 0.3 {GeV/c)"^, implying 
an upper limit of ^ 8% for the increase of the proton charge radius (cf. [|6l|] ). 
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Figure 12. The EMC ratio [Eq. (|)] in '^He (a) and ^^Fe (b) at = iq {GeV/cf. Experimental 
points and the solid lines are as in Fig. 9. The dashed lines are the results of the rescahng formula 
(|5^, calculated adopting in Eq. (|53|) Xa/\n = 1-036 and 1.047 for ^He and ^^Fe, respectively. 



Thus, we assume that {Xa/X^ — 1) is proportional to the nuclear density p{A) in a such 
a way that an increase of 6% is reached only for the heaviest nuclei (namely, at p{A) = 
0.15 /m-3); this corresponds to Xa/Xn = 1.036 and 1.047 for '^He [p{^He) = 0.089 fm-^] 
and ^^Fe [p{^^Fe) = 0.117 fm^^], respectively. In case of the deuteron we assume no swelling 
(i.e., Xd/Xn = 1). The results of the calculations, adopting for the mass scale p^ in Eq. ( ^3)) 
the value 0.6 {GeV/cy as in |^9|, are reported in Fig. 12 for "^He and ^^Fe and in Fig. 13 
for the mirror A = 3 nuclei. It can be seen that the Q^-rescaling approach [Eq. (|52|)] provides 
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a better description of the EMC data at large x for both ^ife and ^^Fe than the convolution 
formula (^31). For ^He and the predictions of the rescaling approach (corresponding to 
XsHe/^p = ^sfj/Xn = 1-020 and XsH^/Xn = Xsff/Xp = 1.027, respectively) provide a possible 
mechanism to achieve a ^ 1.5% deviation of the super-ratio S71emc{x,Q'^) from 1 already 
around x ~ 0.7 0.8. 
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Figure 13. (a) The EMC ratio [Eq. (||)] in ^He (thin lines) and in (thick lines) vs. x at 
= 10 (GeV/c)'^. (b) The super-ratio [Eq. (§)] vs. x. The solid and dashed lines correspond 
to the convolution (^3|) and rescaling ( |52[ ) formula, respectively. In the latter case the values 
X'iRe/Xp = Xsfj/Xn = 1.020 and XzHe/Xn = Xzu/Xp = 1.027 are adopted in Eq. (^). 



4.3 Six-Quark Clusters 

Another mechanism proposed for the explanation of the EMC effect is the formation of 
clusters of six (or more) quarks when two (or more) nucleons are overlapping in the nucleus 
||51| , |52| . This mechanism also provides a softening of the quark distribution in nuclei, since 
the phase space available in a cluster of six (or more) quarks is clearly larger than in a nucleon. 
In what follows we limit ourselves to the case of 6g clusters and we adopt the procedure of 



Ref. ^2[ in order to evaluate the 6g cluster contribution to the nuclear response. 

The main point is to take into account the decomposition (^) of the nucleon Spectral 
Function into a ground {Pq) and a correlated {Pi) parts. Indeed, since two nucleons can 
overlap only in the correlated part P/^, the modification of the convolution formula (^Sf) due 
to the possible presence of 6g clusters can be written as 



N=l''^ ^ N=l "^l ''^ ^ 
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PeX^"'\x,Q% (54) 
where, following Eqs. (|31| ) and (0), one has 

POO POO ]Vf 

fl'iz) = 2nMC^ dE dppP!'{p,E)^==^, (55) 

with i = 0,1. In Eq. (ISl) P^q is the probability to have a six-quark cluster in the nucleus, 
is the normalization of the correlated part of the nucleon Spectral Function, viz. 

/oo POO 
dE dpp^P^{p,E), (56) 



while f/*-^'^''(x, Q^) is given by 



F^''\X, Q^) = 4 E Ir d^^rn Z^m F {z,^) F^ {-^ , Q') , (57) 



where /? = {v?d'^,v?d^ ,u^d'^) = {[nn], [np], [pp]) identifies the type of 6q cluster, f^{zcm) is 
the light-cone momentum distribution describing the center-of-mass motion of the 6g cluster 
in the nuclear medium and -p2^(^,<5^) is the structure function of the 6g cluster. Following 
Ref. p2[, we adopt for f^{zcm) the momentum distribution of the center-of-mass motion of 



a correlated A^A^ pair (with the same quark content of the 6q cluster) as resulting from the 
Spectral Function model of Ref. . In this way we take into account that the 6q bag may 



be not at rest in the nucleus. 

As for -^2^(^,(5^) we follow a simple Q^-independent parameterization proposed in Ref. 



5^ (a) and inspired by quark counting rules, viz. 

FiiO = l« (1 - 0'' + (Es' 0'' \ , (58) 



where the coefficients a and b can be found in Ref. (a). Note that the charge factor c|) 
is different for the various types of 6q clusters, namely: J2j ^"j = 4/3, 5/3, 2 for [nn], [np], [pp] 
clusters. Therefore, the 6q cluster contribution is different in ^He and in ^H, because [nn] 
and [pp] bags have at least different quark content. Note that an additional difference may 
come from different x-distributions in [pn] and [nn], [pp] bags. 

The only remaining parameter is the 6q bag probability P^q in the nucleus. Since the 
probability for two nucleon overlapping is proportional to the nuclear density, we assume Pgg 
to be proportional to the density p{A). We fix the constant of proportionality by requiring 
the best reproduction of the EMC data of ^^Fe, obtaining in this way Pgg ~ 15% in iron 
and Pgq ^ 11% in ^He. We assume no 6q bag in the deuteron, while for mirror A = 3 nuclei 
we get P\pp](inn]) = 6.4% and P^np] = 8.6%. The results of our calculations are reported in 
Figs. 14 and 15. It can be seen that the presence of 6q bags can have an important impact 
on the possible difference of the EMC effect in ^He and ^H, leading to a deviation of the 
super-ratio (0) of ~ 3% already around x ~ 0.7 -h 0.8. 
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Figure 14. The EMC ratio [Eq. (||)] in ^He (a) and ^^Fe (b) at = iq {GeV/cf. Experimental 
points and the sohd hnes are as in Fig. 9. The dotted hnes are the results of the Qq bag formula 
(|5^) , calculated adopting Pgg = 11% and 15% for '^He and ^^Fe, respectively. 



4.4 Color Screening Model 



In inclusive y4(e, e')X reactions the most significant EMC effect is observed at a; ~ 0.5 -^ 0.6. 
This range of x corresponds to high momentum components of the quark distribution in 
the nucleon and therefore the EMC effect is expected to be mostly sensitive to nucleon 
wave function configurations where three quarks are hkely to be close together [l^, . We 
refer to such small size configurations of quarks as point-like configurations {PLC). It is 
then assumed that for large x the dominant contribution to F^lxjQ"^) is given by PLC of 
partons which, due to color screening, weakly interact with the other nucleons. Note that 
due to pQCD evolution F^^x.Q"^) at x > 0.6, Q"^ > 10 {CeV/cY, is determined by the 
non-perturbative nucleon wave function at x > 0.7. Thus it is actually assumed that in the 
non-perturbative nucleon wave function point-like configurations dominate at x > 0.7. 

The suppression of PLC in a bound nucleon is assumed to be the main source of the 
EMC effect in inclusive DIS ||16|, O- Note that this suppression does not lead to a noticeable 



change in the average characteristics of the nucleon in nuclei |T^. To calculate the change of 
the probability of a PLC in a bound nucleon, one can use a perturbation series over a small 
parameter, k, which controls the corrections to the description of the nucleus as a system of 
undeformed nucleons. This parameter is taken to be the ratio of the characteristic energies 
for nucleons and nuclei: 



K 



(Ua) 



AE, 



1 

To' 



(59) 



where (Ua) is the average potential energy per nucleon, {Ua) U>i~ —40 MeV, and AEa 
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Figure 15. (a) The EMC ratio [Eq. (|^)] in ^He (thin hnes) and in (thick hnes) vs. x at 
= 10 (GeV/c)'^. (b) The super-ratio [Eq. (^)] vs. x. The sohd and dotted hnes correspond to the 
convolution (^) and 6q bag (|5^ formula, respectively. In the latter case the values P[pp]([nn]) = 6.4% 
and P[np] = 8.6% are adopted. 



M* — M 0.6 ^ 1 GeV is the typical energy for nucleon excitations within the nucleus. 

The task now is to calculate the deformation of the quark wave function in the bound 
nucleon due to suppression of the probability of PLC in a bound nucleon and then to account 
for it in the calculation of Q^). To this end we consider a model, in which the interac- 

tion between nucleons is described by a Schrodinger equation with a potential V{Rij,yi,yj) 
which depends both on the internucleon distances (besides nucleon spin and isospin) and the 
inner variables i/i and yj, where yi characterizes the quark-gluon configuration in the i-th 



nucleon [|16], |T^, 

In the non-relativistic theory of the nucleus the internucleon interaction V{Rij,yi,yj) is 
averaged over all yi and yj, and the Schrodinger equation is solved for the non-relativistic 
potential U{Rij), which is related to V{Rij,yi,yj) as follows: 



UiR, 



(0iv(?/i)0jv(?/j) I v{Rij,yi,yj,yi,yj) I 0iv(i/i)0iv(z/i)), 



(60) 



yi,yj,yi,yj 



where (pNiVi) is the free nucleon wave function. Using for the unperturbed nuclear wave 
function the solution of the Schrodinger equation with U{Rij), one can treat {U — V)/{Ei — 
Efyf), as a small parameter to calculate the dependence of the probability to find a nucleon 
in a PLC on the momentum of the nucleon inside the nucleus. The quantity Ei introduced 
above is the energy of an intermediate excited nucleon state. Such a calculation allows to 
estimate the suppression of the probability to find a PLC in a bound nucleon as compared to 
the similar probability for a free nucleon. In the DIS cross section the PLC suppression can 
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be represented as a suppression factor 7a (p^) which is multiphcative to the nucleon structure 
function F.^{< x >,Q^) in the LC convolution formula (|37D , viz. [p!6| 

1 1 



(61) 

where Ai?^ = {Ei — E^) ^ M* — M and p is the momentum of the bound nucleon in the 
light cone. 

The X dependence of the suppression effect is based on the assumption that the PLC 
contribution in the nucleon wave function is negligible at x < 0.3, and gives the dominant 



contribution at x > 0.5 [|T6], |31|]. We use a simple linear fit to describe the x dependence 



between these two values of x [^. Using Eq. (plj) for large A at x ~ 0.5 when Fermi motion 



effects are small one can obtains an estimate for Ra in Eq. (|T]) for large A as follows: 

Ra{x) U^o.5~ Mp') ~ 1 + ~ 0.7 - 0.8, (62) 

where (Ua) ~ —40 MeV. Since (Ua) ~ {pa{j)) for A > 12, the model predicts also the 
A dependence of the EMC effect, which is consistent with the data However for the 
lightest nuclei where the Fermi momentum distribution is very steep, the A dependence 
due to the nuclear density is rather oversimplified. The correct estimation requires the 
convolution of Eq. (|6T|) with the structure function of a bound nucleon in Eqs. (|37|) and 



To estimate the suppression factor for large Fermi momenta when the interacting nucleon 
belongs to nucleonic correlations, we use the same formula (|6TD, in which now 7 is defined 
through the virtuality of the interacting nucleon in many-nucleon correlations as follows: 

^ ~ M ■ AEa ' 

\ J J J 

where Mj ^ j ■ M and Mr ^ {j — 1)M are the masses of the j-nucleon correlation and recoil 
(j — l)-nucleon system. 

In Fig. 16 we present the comparison of the predictions of the color screening model 
with the EMC data for ^He and ^^Fe targets. The shaded area is defined by the values 
AE = 0.6 GeV (lower solid curve) and AE = 1 CeV (upper solid curve). 

In Fig. 17 our results obtained for ^He and nuclei assuming the same range of values 
of AE are presented. If one assumes that AE is the same for all nuclei, then even though we 
have large uncertainties for the EMC ratio [see Figs. 17(a) and 17(b)], the nuclear effects 
largely cancel out in the super-ratio (0) [see Fig. 17(c)]. However the assumption that AE 
is independent on the specific nucleus is a clear oversimplification, since the bound nucleon 
excitation in NN correlations does depend on the isospin of the A^A^ pair. Indeed, as it is 
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Figure 16. The EMC ratio in ^He (a) and ^^Fe (b) at = iq (GeV/c)'^. Open dots are data 
from Ref. while in (b) the full squares are from Ref. [^^. The dashed lines are the results 

of Eq. (|37|) calculated adopting the parameterization of the nucleon structure function F2^{x,Q'^) 
of Ref. [Q, including nucleon resonances, and the model of Ref. [33| for the nucleon Spectral 
Function P^{k,E). The shaded area corresponds to the prediction of the color screening model 
with AE = 0.6 GeV - lower solid curve and AE = 1 GeV - upper solid curve. 



shown in Fig. 10, one has more attraction in isosinglet than in the isovector pairs and the 
spatial distributions for different isospin states may substantially differ. Therefore, one can 
expect that AE is smaller in isotriplet states as compared to the isosinglet cases. To estimate 
the upper limit of uncertainties due to the expected isospin dependence of bound nucleon 
excitation, we decompose the A^A^ correlation into an isosinglet and an isotriplet contribu- 
tions. Then we assume for isotriplet states the minimal value of AE = 0.6 GeV, whereas 
for isosinglet states we assume no suppression at all. The prediction of this approximation 
is shown as the dotted curve in Fig. 17(c). 

The isospin dependence of the EMC effect emerges naturally also in the quark-diquark 
model of the nucleon. Within this model one expects a different degree of suppression for 
valence quark and quarks in the diquark. In Fig. 17(d) we have presented the predictions of 
the color screening model within the quark-diquark picture, in which we have considered two 
extreme cases: the first one when the suppression occurs only for the valence quarks (thin 
solid curve) and the second one when only the quarks in the diquark are suppressed (thick 
solid curve). To estimate the extent of the suppression we assume AE ^ M^^ — M ~ 0.3 GeV 
which roughly corresponds to the quark helicity-flip excitation in the nucleon. 

As follows from Figs. 17(c) and 17(d) all these approximations produce at most a 2 -=-3% 
deviation in the super-ratio (|]) around x ~ 0.7 -r- 0.8, which may be considered as an upper 
limit of the uncertainties due to the EMC effect in mirror A = 3 nuclei within the color 
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screening model. 



5 Extraction of the n/p Ratio 



We now use the the results of the previous Sections to address the issue of the extraction of 
the neutron to proton structure function ratio F2 {x , Q"^) / F2 {x , Q"^) from the measurements 
of the ratio of the mirror A = 3 structure functions, defined as 



nF{x,Q^) 



Fl^\x,Q^) 



(64) 



From Eq. one gets 



F^{x,Q^) 



2snEMcix,Q^) -nF{x,Q^' 



and correspondingly the uncertainty on the extracted n/p ratio is given by 
A (F^/Fi) 3 SUEMcix, Q') TZpix, Q^) 



(65) 



F^/F 



[27^^(x, g2) _ suEMcix, g2)] [2snEMc{x, Q') - 7^^(x, q^) 



\ 



A{sn 



EMC) 



sn 



EMC 



+ 



A(^ 

7^F 



(66) 



One can easily see that in the r.h.s. of Eq. (|66|) the quantity in front of the square 
root provides a factor of ~ 4 at x > 0.7; thus, even a small uncertainty on the super-ratio is 
largely amplified in Eq. (^), yielding a non- negligible uncertainty on the extracted n/p ratio 
at large x. Following Ref. the total experimental error in the DIS cross section ratio of 

and ^He is likely to be < 1%. Therefore in Eq. (|66D we assume that A(71e)/'R-f = 1% 
and from the results of the previous Sections we consider that A.{S71emc) / STZemc = 3% 
for X > 0.6, leading to a total uncertainty of ~ 12% in the extracted n/p ratio already at 
X ^ 0.7. In Fig. 18 we have reported the x-dependence of the expected accuracy for the 
extraction of the ratio -F^(x, Q^)/F2 (x, Q^), and the shaded areas include the combination of 



all the effects discussed above. Moreover we have used the CTEQ flOl set of PDF's (lower 
shaded area) and the modified CTEQ parameterization (upper shaded area), obtained from 
the CTEQ one as described in Sec. 3, in order to reproduce a. n/p ratio approaching 3/7 at 
X — > 1. Figure 18 demonstrates that although mirror A = 3 measurements will significantly 
improve the existing accuracy of the neutron structure functions at large x, they may not 
provide 3cr separation for the two predictions of the n/p ratio having limiting values of 1/4 
and 3/7 at x ^ 1. 

In Ref. ||2^, ^ it was suggested that, once the ratio 71f{x,Q'^) is measured, one can 



employ an iterative procedure to extract the n/p ratio which can almost eliminate the effects 
of the dependence of the super-ratio STZEMcix^Q"^) on the large-x behavior of the specific 
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structure function input (see Fig. 6). Namely, after extracting the n/p ratio assuming a 
particular calculation of STZemc, one can use the extracted neutron structure function to 
get a new estimate of STZemc, which can then be employed for a further extraction of the 
n/p ratio. Such a procedure can be iterated until convergence is achieved and self-consistent 
solutions for the extracted F^^/Ff and the super-ratio SIZemc are obtained. In Ref. p6 



the numerical estimate of the iteration procedure was performed within the VNC model 
and a good convergence was achieved for x up to ^ 0.8. However this result depends on the 
assumed validity of the VNC model in the considered range of values of x. To check how 
well the iteration procedure will work in case of other models of the EMC effect, we have 
considered the following two examples. 

First, let us consider in Eq. (|65D the nuclear structure function ratio 71f{x,Q^) which 
results from the use of the modified SLAC parameterization of F<^{x, Q"^) and the inclusion of 
the effects of possible 6g bags within the VNC model. Then we apply the iteration procedure 
assuming for the super-ratio S71emc{x, Q"^) the convolution calculation corresponding to the 
Donnachie-Landshoff {DL) fit |^ of F.^{x,Q'^), which provides a n/p ratio equal to 1/3 as 
X 1. Figure 19 demonstrates that a consistency is achieved between the n/p ratio used as 
the input and the extracted one. However, for the calculation of 71e{x, Q"^) we started from 
the modified SLAC parameterization which goes to 3/7 as a; —> 1. 

In the second example we calculate the nuclear structure function ratio IZf^x^Q"^) us- 
ing the modified CTEQ parameterization within the LC approximation adopting the color 
screening model for the EMC effect (only with valence quark suppression in the quark- 
diquark picture). To do the iteration we start with the super-ratio S71emc{x, Q"^) calculated 
within the LC approximation without EMC effects using the CTEQ parameterization for 
the nucleon structure function. Figure 20 demonstrates that the iteration diverges already 
at values of x (~ 0.7) smaller than the ones obtained in E^, where only the VNC model 
was used. Below x ^ 0.7 the iteration procedure converges to a value of the n/p ratio which 
differs from the "exact" one (used in Re{x,Q'^)) exactly by the amount of the EMC effect 
which is implemented in the calculation of Rf{x, Q^)- 

Both examples illustrate that the iterative procedure cannot improve the accuracy of the 
extraction of the n/p ratio as estimated in Fig. 18. 



6 Conclusions 

In conclusion, we have performed a detailed theoretical investigation of the EMC effect in 
^He and targets. First, we have considered first the convolution approximation within the 
VNC model and the LC approach. The differences in the motion of protons and neutrons in 
mirror A = 3 nuclei, resulting from the spin-flavor dependence of the nuclear force, have been 
taken into account. We have investigated various sources of uncertainties in the estimate 
of the super-ratio of the EMC effects in mirror A = 3 nuclei, like: i) charge symmetry 
breaking terms in the nucleon-nucleon interaction; ii) finite corrections to the Bjorken 
limit; iii) the role of different prescriptions to relate the invariant nucleon Spectral Function 
to the non-relativistic one, required to ensure baryon number conservation; iv) the role 
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of the x-shape of parton distribution functions; and v) the differences between the VNC 
model and the hght cone formahsm. 

Within convolution approach, in which no modification of the bound nucleon is consid- 
ered, the deviation of the super-ratio (§) from unity is predicted to stay within 1% only for 
X < 0.75 in close agreement with Ref. |^ and in overall agreement with Ref. (which by 
the way neglect the effects (i)-(iv)), where 1% deviations were found from an average value 
of the super-ratio equal to ~ 1.01. 

We have further argued that the previous estimate cannot be considered as definitive 
for the purpose of extraction of the F^/Ff ratio, since it is derived using just one of the 
many models of the EMC effect, which has in particular a number of problems in describing 
the nuclear data, namely the underestimation the EMC effect at x > 0.6. Hence we have 
provided a detailed analysis of the super-ratio within a broad range of models of the EMC 
effect, which take into account possible modifications of the bound nucleons in nuclei, like: 
i) a change in the quark confinement size (including swelling); ii) the possible presence 
of clusters of six quarks; and iii) the suppression of point-like configurations due to color 
screening. Our main result is that one cannot exclude the possibility that the cancellation 
of the nuclear effects in the super-ratio may occur only at a level of ~ 3%, resulting in a 
significant uncertainty (up to ~ 12% for x ~ 0.7 ^ 0.8) in the extraction of the free n/p 
ratio from the ratio of the measurements of the ^He and DIS structure functions. Such 
an uncertainty is comparable to the ~ 18% difference between the n/p predictions having 
limiting values of 3/7 and 1/4 at x 1, which characterize the pQCD and Feynman models. 
Another consequence of the use of a broad range of models for the EMC effect is that the 
iteration procedure cannot in general improve the accuracy of the extracted n/p ratio. 

It is however important to note that despite of such restrictions the mirror A = ?> measure- 
ments will provide an unprecedented accuracy in the extraction of the neutron DIS structure 
function. Thus such measurements are very much welcomed. It is however very important 
to complement them with the measurements of semi-inclusive processes off the deuteron, in 
which the momentum of the struck nucleon is tagged by detecting the recoiling one. Impos- 
ing the kinematical conditions that the detected momentum is low {p < 150 MeV/c), which 
means that the nucleons in the deuteron are initially far apart [p!7| , p3| , p4| , it is possible to 
minimize significantly the nuclear effects. Furthermore, all the unwanted nuclear effects can 
be isolated by using the same reaction for the extraction of the proton structure function 
by detecting slow recoiling neutrons and comparing the results with existing hydrogen data, 
as well as by performing tighter cuts on the momentum of the spectator proton and then 
extrapolating to the neutron pole. 
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Figure 17. The EMC ratio in ^He (a) and in (b) vs. x at = 10 (GeV/c)'^. The super- 
ratio [Eq. ^]) vs. X (c, d). In (a)-(c) the dashed curves correspond to the LC formahsm ( ^7\ ) 
using for the nucleon structure function the parameterization of Ref. |^], which includes nucleon 
resonances. In (a)-(c) the shaded areas correspond to the predictions of the color screening model 
for isospin independent values of AE, namely for AE = 0.6 GeV - lower (upper for (c)) solid curves 
and AE = 1 GeV - upper (lower for (c)) solid curves. In (c) the dotted line corresponds to the 
prediction of the isospin dependent screening model described in the text. In (d) the thin and thick 
solid curves correspond to the predictions of the screening model within the quark-diquark picture, 
when only the valence quarks and the quarks in the diquark are suppressed, respectively. 
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Figure 18. The expected accuracy for the extraction of the neutron to proton ratio 
^2^(3;, Q^) vs. X at = 10 (GeV/c)'^. The lower and upper shaded areas correspond 

to the CTEQ and modified CTEQ parameterization described in the text. The full dots are the 
NMC data points as given in Ref. |45]. 



X 



O 



0.7 



O 0.6 



^ 0.5 
0.4 
0.3 

CVI 




. 2 LI I I I I I I I I I I I I I I I I I I I I I I I I I I L 



0.3 



0.5 



0.7 



0.9 



Figure 19. The neutron to proton ratio F^{x,Q'^) / F^{x,Q^) vs. x aX Q"^ = 10 {GeV/cf. The 
full dots are the NMC data points as given in Ref. EH]. The solid line represents the DL 



parameterization [33| of the n/p ratio. The open squares are the results of the extraction of the 
n/p ratio adopting the the convolution formula (23) as described in the text. The error bars are 
calculated via Eq. (B^). 
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Figure 20. The x dependence of neutron to proton structure function ratios-F2n/-^2p- The thick 
sohd curve corresponds to the n/p ratio used to calculate Rp- The thin solid curve is the ratio used 
as input for the iteration procedure. The dashed curves correspond to the extracted ratio obtained 
after subsequent iterations. 
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